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Abstract. We consider a simple model of communities interacting via bilinear terms. After 
analyzing the thermal equilibrium case, which can be described by an Hamiltonian, we introduce 
the dynamics that, for Ising-like variables, reduces to a Glauber-like dynamics. We analyze and 
compare four different versions of the dynamics: flow (differential equations), map (discrete¬ 
time dynamics), local-time update flow, and local-time update map. The presence of only 
bilinear interactions prevent the flow cases to develop any dynamical instability, the system 
converging always to the thermal equilibrium. The situation is different for the map when 
unfriendly couplings are involved, where period-two oscillations arise. In the case of the map 
with local-time updates, oscillations of any period and chaos can arise as a consequence of the 
reciprocal “tension” accumulated among the communities during their sleeping time interval. 
The resulting chaos can be of two kinds: true chaos characterized by positive Lyapunov exponent 
and bifurcation cascades, or marginal chaos characterized by zero Lyapunov exponent and 
critical continuous regions. 


1. Introduction 

The fact that the speed of physical interactions is finite, as for light or sound, has well known 
crucial consequences in all physics, and it is the ultimate reason for the existence of waves. For 
example, if we artificially assume the speed of light as infinite in the Maxwell equations, the 
magnetic field would be zero and we would not have electromagnetic waves. 

The situation is quite different when the equations contain non linear terms, which can result 
in erratic or even chaotic trajectories. Notice that, in the absence of mechanical kinetic terms, 
which is our framework, the presence of non linear terms is exclusively due to the presence of 
non bilinear interactions among the components of the system. For example, in the Kuramoto 
model M-m, sinusoidal interactions lead to chaos m-m if the coupling constant K is larger 
than some threshold Ke¬ 
lt should be emphasized that the presence of non linear terms is not always sufficient for the 
emergence of chaos. However, it has been observed that the introduction of time-delays among 
the components of a system augments the chance that the system develops chaos [9]-[13]. The 
idea at the base of such a mechanism is that time-delays among the system components makes 
it difficult for the system to set in a steady state. In particular, for linear flows, a fixed delay 
among the system components corresponds to a finite speed of interactions, which in turn leads 
to oscillations or waves. For non linear flows a delay can have more dramatic consequences, with 
chaos scenarios more probable. 

In |17j . we have introduced a simple discrete-time two components dynamical model, 
characterized by local-time updates. Local time updates means that the each component updates 


its status only at certain times t. In the specific model analyzed in HZl, a component updates 
(sleeps) its status only at odd times t = 1,3,..., (even times t = 0,2,...), while the other 
component updates (sleeps) its status only at even times t = 0, 2,... (odd times t = 1,3,...). 
Remarkably, as we have shown in ini, for certain coupling values, this simple system exhibits 
chaos even if it is characterized only by bilinear interactions. In fact, when the couplings are in 
the steady state regime, the system can be described via a Hamiltonian which is quadratic in 
the system components. We are not aware of a similar case in the literature. Our understanding 
here is that there are three factors generates chaos: i) the presence of unfriendly couplings, 
ii) the discrete-time nature of the dynamics, and Hi) the existence of at least two interacting 
communities which rearrange their configurations at alternate times, i.e. via local-time updates. 
Notice that (in) is not just a delay; it also implies that each component has a sleeping period. 
It is during such a sleeping period that the component accumulates a tension with respect to 
the other component when the reciprocal coupling is unfriendly. We are in fact dealing with a 
sort of dynamical frustration. 

In this paper we review some aspects of the model which were only partially treated in m- 
In particular, we compare the analogous flow model with the present map and show that the 
above conditions (i)-(iii) are in fact necessary, and calculate the maximal Lyapunov exponent. 

The reader not interested in the derivation of the dynamic equations may skip the next 
Section. 


2. The model in the flow, map and local-time versions 

Here we consider a two components system. Generalization to an arbitrary number of 
components will be reported elsewhere. We consider two communities [TT] of agents and 
with cardinalities and The agent with index i can be in two possible status, 

fji = ±1, with i G UA/'^^^ According to the sign of the couplings, friendly or unfriendly, 
each agent i tends to follow or anti-follow its neighbors, by minimizing or maximizing the term 
cJi^^cTj, where runs over the set of neighbors of L It is necessary to distinguish between 
intra- and inter-couplings, therefore we introduce the 2x2 matrix jLd) j(2:2) 

the intra-couplings, and the inter-coupling. In the most general formulation we 

should introduce also the 2x2 matrix defined as the set of coupled spins (i,j) within the 

same community (intra), or between the two communities (inter). Finally, we introduce a global 
factor f3 that rescales all the couplings in 

Let us indicate the time by t, which is supposed to range in the set of natural numbers N for 
the map, while it can take any non negative real number for the flow. Inspired by the Glauber 
dynamics |18j . for flow and map we introduce the following transition rate probabilities for the 
spin with state ai to jump to the state 


w 


(FLOW/MAP) 


(o-j cr';t) = < 




Let (E,0) be any periodic partition of N, i.e., EuO = N, EnO = 0, and E and O contain 
both infinite elements of N. For instance, E and O can be the set of even and odd numbers, 
respectively, and we will assume this choice in all the next examples. We now introduce the 
following local-time transition rate probabilities for the spin with state CTj to jump to the state a[ 





( 2 ) 
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if mod[i,2]/0, 

l+a' tanh(/3j(i'i) E^,(i,^.)gr(i. 2 ) ^j) 
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for i G N'^^\ if mod[f, 2] = 0, 

i G if mod[i, 2] / 0, 

1+<t' tanh(/3j(2,2) E^^(,,^^)gr(2.2) ^,-+/3J'2,i) E^,(,,^.)gr(2.i) <^i) 
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(LOCAL-TIME MAP) ^ ^ ^ 


for i G t G E, 

for t G E, 

L 1 2 

l+f7'tanh(/3j(2.2) E^,(,_^.)gp{2,2) -Tj+/3J(2.i) E^,(,_ .)gp{2,i) 

2 ' 

for i G t € O. 


(3) 


From the viewpoint of modeling, Eqs. 0, (§ , or (|^, are justified as they make each spin to 
follow the majority of its intra- and inter-neighbors and, thanks to the presence of the functions 
tanh(-), the rates are non negative and normalized at any time ^'it) = l/(Time 

Unit). From a deeper viewpoint, Eqs. ([^, Q, or (|^, are based on the fact that, as we will 
see soon, in the case of positive couplings they lead to Boltzmann equilibrium governed by only 
quadratic interactions. More precisely, in the case of positive couplings these forms guarantee 
that at equilibrium the system satisfies the principle of detailed balance and the principle of 
maximal entropy for any quadratic interactions. 

We formalize the discrete-time probabilistic dynamics induced by Eqs. 0, 0 , or 0 , or 
as follows. Let N = + N^'^\ Let us introduce the spin vector cr = (cti, ... ,(7^), and the 

associated probability vector p(<t; t), i.e., the probability that the system is in the configuration 
a at time f G N. The master equation for flow and map (local-time or not) reads 


clp(FLOW) 

d{at) 


/ 

O' 


/ 

O' 


p(MAP) (cr; t -|- 1) — p('^A.P) 
a 

+ J2p^^^^Ho-';t)W{a' cr), 

/ 

O' 

where we have introduced the global transition rates 

W{(T-^a')= 0 wiai^a-), 

ieAr(i)uAr(2) 


/ 

O' 
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( 6 ) 










where w(ai —)• can be of the type Eqs. 0.0 , or Q , and where a/2 > 0 may be interpreted 

as the rate at which, due to the interaction with an environment, a free spin (= 0) makes 
transitions from either state to the other. As we have proved in m, it is necessary to impose 
the bound a < 1 for p{(T]t) to be a probability at any time t. By using Eqs. (§-® it is easy 
to check that the stationary solutions p(cr) of Eq. ([^ are given by the Boltzmann distribution 
p{cr) oc exp[—f3H{cr)], where 

H= — ^ ataj — ^ UiUj — ^ (JiOj- (7) 

(i,i)er(i.i) (j,i)er(2.2) (jj)gr(i.2) 


Of course, the existence of a stationary solution p(cr) does not represent a sufficient condition 
for equilibrium. In fact, asymptotically the system can reach non-point-like attractors, and even 
aperiodic or chaotic regimes. However, as we have discussed in the Introduction, the fact that 
for certain values of couplings the system can be described via the Hamiltonian 0 , means that 
we are dealing with only bilinear interactions. Here, the non linearities that lead to chaos are 
caused by the local time updates. 

From now on, we will omit the s uffix es (FLOW), (MAP) and (LOCAL-TIME). It will be clear 
from the context what we are referring to. Eqs. 0-0 define the microscopic dynamics from 
which one can derive the macroscopic dynamics, i.e., the dynamics for the order parameters 


= '^p{a-,t)^ ai. (9) 


ieA/'(2) 


3. Mean Field Limit 

The mean-field limit is defined by the settings = (^ 2 *)’ = (^ 2 *)’ = 

and the replacements —>• —)• —)• 

A^0))/(2A'(^^A^(^)). We parametrize the size of the two communities as 


Ar(i) = Np^^\ Ar(2) = iVp(2)^ pW + p(2) = 


Let us introduce the matrix J 


J = 


7(1,1) J(1’2) 

2p(l) 

j(2,l) j(2,2) 

2p(2) 


( 10 ) 


( 11 ) 


As shown in El, for —7- 00 we obtain the following deterministic evolution Eqs. for the order 
parameters 0-0 


= tanh (/3j(i’i)x(i)(t) + /3J(i’2)a;(2)(t)) 
= tanh (/3j(2-2)a;(2)(i) + ^j(2,i)^(i)(i)^ 



(12) 


(13) 


3.1. FLOW 






3.2. MAP 


^ (f(j(2,2)a:(2)(t) + ,3j(2.i)^(l)(t)) _ ^(2)(t)_ (15) 


3.3. LOCAL-TIME FLOW 


dx^^\t) 

d{at) 


0 , if mod[f, 2 ] / 0 , 


tanh _|_ ^j(i 

2 )x( 2 )(t)) 


dx^‘^\t) 

d{at) 


0 , if mod[f, 2 ] / 0 , 
tanh 


(16) 

x^^\t), if mod[f, 2 ] = 0 , 


(17) 

x^‘^\t), if mod[f, 2 ] = 0 . 


3 . 4 . LOCAL-TIME MAP 

x^^\t + 1 ) — 
a 

x^‘^\t + 1 ) — x^‘^\t) 
a 


' 0 , t £0, 

^ / ~ ~ X 

tanh — x^^\t), f G E, 

' 0 , f G E, 

^ ~ ~ X (^9) 

tanh —x^'^\t), t G O. 


3.5. Stationary solutions 

When the couplings are positive, there is little difference among the four dynamics, and 
they all asymptotically reach equilibrium according to the Boltzmann distribution p(cr) oc 
exp[—/377iiif(cr)], where is the mean-field analogous of Eq. Q, and x^^\t) and x^'^\t) 

tend, for t —)■ 00 , to the stationary solutions of Eqs. i.e.. 


x^^'^ = tanh -|- , 

= tanh -|- . 


( 20 ) 


Eqs. (20) represent a particular case of the general result derived in m valid for n interacting 
communities at equilibrium. In particular, one can check that Eqs. (20) give rise to second order 


phase transitions whose critical surface is determined by the condition 


det (1 — 


/3j) 


= 0 . 


( 21 ) 


In general, (x^-^^ x^^^) = (0,0) is stable when the eigenvalues of j3J are inside the interval (—1,1) 
(disordered phase, or no consensus), otherwise the system reaches a spontaneous ordered phase 
(x^^^x*-^^) 7 ^ (0,0) (frozen phase, or consensus). However, Eq. (21) is exact only for the flow 
and the map. In general, due to the sleeping time, the critical value of /3 can differ from the 
one given by Eq. (HI)- This was erroneously not stressed in m- The correct way to to derive 
the exact critical value of j3 must use the full expression of the Jacobian, which is the subject of 
Section 5. 















4. Comparing flow, map and local-time versions 

From now on we shall always assume a = 1 and 

In this Section we compare the four dynamics numerically. 

In Figs (la)-(2b) we compare the normal flow and map. In Figs, (la) and (lb) we consider a 
situation with positive couplings at small and large values of /3; in Figs. (2a) and (2b) we consider 
a situation with positive and negative couplings at small and large values of (3, respectively. It 
is evident that, if the couplings are positive, the flow and the map remain roughly close at any 
time, and tend to coincide for t —>■ oo to the same stationary solution provided by Eq. (20). 
When instead some of the couplings are negative, at small /3 the map tends to make damping 
oscillations around the flow (which for the present case is null), whereas, at large /3, the map 
tends to make regular oscillations around the flow solution (which for the present case is null). 


(a) 

P=0.3, 




Figure 1: (Color online) Plot of and solution of the flow, Eqs. (12)-(13), and the 

map, Eqs. (14)-(15) (here with a = 1), for positive couplings; below (Panel a) and above (Panel 
b) the threshold j3c- Initial conditions are the same in all cases. The lines of the map are guides 
for the eyes. 
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(a) 

M.3, J"''W<“>=-1, 
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(b) 



Figure 2: (Color online) Plot of and x solution of the flow, Eqs. ((l^-((l^, and the 

map, Eqs. (@-([151) (here with 0 = 1), for a case with positive and negative couplings; below 
(Panel a) and above (Panel b) the threshold Pc- Initial conditions are the same in all cases. The 
lines of the map are guides for the eyes. 


In Figs. (3a)-(3c) we compare the flow and map with local-time updates for a case with 
positive and negative couplings and three values of /3, respectively. Fig. (3a) provides the 
behavior for /3 = 1.5 where the map follows damping oscillations and both the map and the flow 
tend to the same stationary value provided by Eq. (20). Eig. (3b) provides the behavior for a 
relatively larger value of /3, /3 = 2.5, where the map follows regular oscillations of period 2, while 
the flow tends to a stationary value. Fig. (3c) provides the behavior for /3 = 4.5, where the map 
follows regular oscillations of period 4, while the flow tends to a stationary value. Increasing 
further /? leads to larger and larger periods (not shown). 
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Local-Time (b) 

[5=2.S,j“-‘’=-l.j“=0.4,j'''">=j‘^'‘>=0.8 
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Local-Time (c) 
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Figure 3; (Color online) Plot of and solution of the local-time flow, Eqs. (16)-(17), 

and the local-time map, Eqs. 


(here with a = 1), (3 = 1.5 (a), /3 = 2.5 (b), and /3 = 4.5 
Initial conditions are the same in all cases. The lines of the map are guides for the eyes. 
















































































































Phase Diagram Local-Time MAP (a) 
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Phase Diagram Local-Time MAP (b) 
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Figure 4: (Color online) Plot of (/3) (red dots), solution of the local-time map, Eqs. (18)-(19) 
(here with a = 1) as a function of /3, with T high enough to remove temporary transients. In 
this case T = 10^. We plot also the maximum Lyapunov exponent A(/3) (black highly irregular 
line) calculated from Eqs. (25)-(27). Each value of /3 is in correspondence with a random initial 
condition for a total of 10^ samples. Panel (a) is a case with bifurcation cascade. Panel (b) is a 
case with marginal chaos. 


5. Chaos 

The case analyzed in Figs. (3) shows that, for large values of /I, the local-time map can lead to 
erratic and, possibly, chaotic trajectories. In general, for each dynamics, the major quantity of 
interest is the Jacobian where l,m£ {1,2}: 


5.1. FLOW 


Ml,m) 
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( 22 ) 



























5.2. MAP 


'Jt+1 — 



5.3. LOCAL-TIME FLOW 


Ml,m) 


di^rn, mod[t + Z - 1,2] / 0, 


1 — [x 




mod[i + /-l,2] =0, 


5 . 4 . LOCAL-TIME MAP 


Ml,m) 

'-'Z+1 


di^rn, mod[i + Z - 1,2] / 0, 


1 — (x 


„(0 


/3j(z,m)^ mod[t + /-l,2] =0. 


(23) 


(24) 


(25) 


The knowledge of allows to calculate the maximum Lyapunov exponent A. It turns out 

that in the map case (local-time or not), due to the fact that - via recurrence - we have the 
explicit solution of the Eqs., the evaluation of A is easier if compared to the flow case. Since the 
flow does not show any signs of chaos or erratic trajectories, we can limit the evaluation of the 
Lyapunov exponent to the map, where we can use 

^ (26) 

At 

where At is the largest eigenvalue of the following matrix 


t 

Qt = LtL\, where Lt=Y\Jt' (27) 

t'=0 

Given a set of couplings, we let the system to evolve toward high enough values of t = T in 
order to remove temporary transients and repeat the numerical experiment for several values of 
/3, each /3 being associated to a random initial condition. Hence, we plot {xP\T),x 0)(T)) as 
functions of /3, and we indicate these functions as {xip\/3), x'^\/3)). In our examples and for the 
range of [3 considered, T > 10^ turns out to be enough high. In fact, the variable /3 plays the role 
of a time-scaling: the higher /3, the higher the necessary transient T. In general, the functions 
(x^^(/3), Xy^(/3)) look multi-valued functions due to the existence of bifurcation points. There 
exist two kinds of bifurcations: doubling period and phase transition. However, it is clear that 
in either case, doubling period or phase transition, the presence of bifurcations increases the 
chance to develop chaos. Therefore, for our aims here, it is more interesting, as well as highly 
more efficient, to show the plots that include all the bifurcation points. Figs. (4a)-(4b) show 
two different scenarios. The case of Panel (a) is a “classical” bifurcation cascade [I], where, 
after a threshold flchaos, bifurcation takes place increasingly, up to windows of stability. In the 
case of Panel (b) we observe a different situation which we call “marginal chaos”. In this case, 
after a threshold (3 Marginal, the dynamics is characterized by totally erratic trajectories with 
no period, up to windows of stability. Strictly speaking, even out of the windows of stability, 
the system is not chaotic, since the maximum Lyapunov exponent is A = 0. However, in such 
regions, the system turns out to be critical over continuous intervals; a situation which is not 
typical in models characterized by discrete symmetries like the Ising model. 










6. Conclusions 

We have analyzed and compared four kinds of dynamics of a simple two component Ising-like 
system: flow, map, local-time flow, and local-time map. Remarkably, even if only bilinear 
interactions are present, the local-time map gives rise to erratic trajectories and, depending on 
the set of couplings, two chaotic scenarios take place: bifurcation cascades or “marginal chaos”, 
i.e., criticality extended over continuous intervals of (3, the time-scale parameter. The analogous 
local-time flow does not present any of such behaviors, not even oscillations, confirming all the 
scenarios we have already discussed in m- We stress that this local-time update is not just 
a time-delay: here we have also a sleeping time during which, in the case of some unfriendly 
couplings, there is an accumulation of frustration among the system components, giving rise to 
instability and chaos, even for a system characterized by only bilinear interactions. As we have 
discussed in m, we believe that maps with local-time updates are a quite common feature in 
complex systems. Urgent issues will be to understand how this scenario generalizes to large 
systems with many, say q components, how the probability to have chaos changes with q, and if 
the mean-held (or more precisely fully connected) picture remains robust. 
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